Anti-lensing: the bright side of voids 
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More than half of the volume of our Universe is occupied by cosmic voids. The lensing magnifica- 
tion effect from those under-dense regions is generally thought to give a small dimming contribution: 
objects on the far side of a void are supposed to be observed as slightly smaller than if the void were 
not there, which together with conservation of surface brightness implies net reduction in photons 
received. This is predicted by the usual weak lensing integral of the density contrast along the 
line of sight. We show that this standard effect is swamped by a relativistic Doppler term that is 
typically neglected. Contrary to the usual expectation, objects on the far side of a void are brighter 
than they would be otherwise. Thus the local dynamics of matter in and near the void is crucial 
and is only captured by the full relativistic lensing convergence. There are also significant nonlinear 
corrections to the relativistic linear theory, which we show actually under-predicts the effect. We 
use exact solutions to estimate that these can be more than 20% for deep voids. This remains an 
important source of systematic errors for weak lensing density reconstruction in galaxy surveys and 
for supernovae observations, and may be the cause of the reported extra scatter of field supernovae 
located on the edge of voids compared to those in clusters. 



Introduction Lensing phenomena are measured not 
only around virialized clusters of galaxies but also 
through and around unvirialized cosmic voids, which oc- 
cupy well above half the volume of the Universe. Here we 
show how the standard lensing magnification effect can 
be overwhelmed by relativistic corrections to the size and 
brightness of sources in and near voids. 

Magnification in the linear approximation The 
lensing magnification effect can be expressed in terms of 
the convergence k, which corrects the background angu- 
lar diameter distance (d A ) 



d A (z) = d A (z)[l - k(z)]. 



(1) 



The convergence in a perturbed ACDM universe is usu- 
ally given as a line of sight integral over the density con- 
trast 5, 



ft = Kg 



d X (XS X) x(l + z)5( X ), (2) 
Xs 



where d^ = dz/H — —dm x is the comoving distance, r\ 
conformal time and S denotes the source. In fact the full 
relativistic expression is p] E] 



ft = fty2$ + K v + Ksw + K/, 



(3) 



where the Sachs- Wolfe term ksw is given by the dif- 
ference in gravitational potential $ between source and 
observer, and ft/ is a line of sight integral over $ and 
its conformal time derivatives <&',<!>". These two terms 
are sub-dominant [T], and we will not discuss their de- 
tailed form, although we do include them in our nu- 
merical calculations below. (The perturbed metric is 
ds 2 = a 2 [-(1 + 2$)d?/ 2 + (1 - 2$)da: 2 ] .) 



The usual form 
term on the right of 
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is an approximation to the first 
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The screen-space Laplacian is Vj_ = V 2 — (n • V) 2 — 
2x -1 n • V, where n is the unit direction from the source. 
The radial derivatives lead to terms proportional to $, 
$' and [2], which are much smaller than the term 
V 2< £ on the sub-Hubble scales of interest. Thus in Q, 
we may replace V 2 ^ by V 2< £, which is given in terms 
of the density contrast 8 by the Poisson equation. The 
general relativistic Poisson equation involves also the pe- 
culiar velocity (vi = div): 



V 2 $ 
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(5 - 3aHv) 



(& + aH<£). 



(5) 
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By ([6| , aHv is of order <I> and may be neglected in |5]) on 
the relevant scales. Then Q reduces to the usual lensing 
term pj). For an under-density, S < 0, so that Kg < 
if the underdensity is the dominant structure along the 
line of sight. Then ([I]) implies that the angular distance 
should be larger than the background value for a fixed 
z, and objects should consequently be observed to be 
smaller. Since surface brightness is conserved in lensing, 
the total number of photons arriving from the object per 
unit time should be fewer. We can quantify this using the 
change to the distance modulus, Am = 51og 10 d A /dA, so 
that a positive Am corresponds to a fainter source. 
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FIG. 1: Change in distance modulus due to a 50 Mpc radius void, vs observed redshift, based on: relativistic convergence 
Q (solid), usual weak lensing formula |2| (dotted), Doppler term Q (dot-dashed), exact model (dashed) (top panels). There 
is a brightening for objects on the far side of the void, whereas the usual weak lensing predicts a dimming of much smaller 
magnitude (see insets). The left panel shows the effect for a small void well in the linear regime, and the right panel shows a 
very deep void where non- linear contributions increase the effect. The bottom panels show the void density contrasts. (Asterisk 
marks the far edge of the void.) 



The usual formula is a good approximation to three of 
the terms in Kg s» Ky 2 * ^ K sw, K i- The remaining 
Doppler term, which arises from a shift in the redshift 
from its background value, 



1 - 



XsH s 



v s ■ n, 



(7) 



is typically ignored - but it cannot be neglected as em- 
phasised by pQ , and as seen in Fig. [I] 

Figure [I] (left) shows the correction to the distance 
modulus for a spherical void of radius 50 Mpc in the 
linear regime, <5min — — 0.05, located at z — 0.1. It is 
clear that k$ predicts a completely wrong magnitude for 
sources in or near the void - and furthermore it predicts 
the wrong sign. By contrast, k v gives a very good ap- 
proximation to the full relativistic k. Near the far edge of 
the void there is a significant positive magnification sig- 
nal: objects are brighter than they would be otherwise, 
an effect which extends far beyond the edge of the void. 
This is the opposite effect one expects based on a naive 
prediction using the usual lensing formula. We also show 
the prediction of an exact model of the void, using the 
Lemaitre-Tolman-Bondi (LTB) solution (see below). It 
demonstrates that the linear relativistic k is accurate for 
this amplitude of void. Since there is no background for 
the LTB case, the effect is not due to peculiar velocity, 
but rather to the extra redshifting of photons as they 
pass through a region of higher expansion rate and non- 



zero shear. The effect is stronger for a deep void with 
<5min = —0.95 (Fig. [I] right panel), where the linear ap- 
proximation underestimates the anti-lensing effect - we 
consider this in more detail below. 

Modelling voids via nonlinear solutions Real 
voids typically have S < —0.8 [3], and we need to ex- 
tend the relativistic perturbative analysis to deal with 
such voids. We can gain some insight via exact solutions 
of the Einstein field equations, where a void region is em- 
bedded in a homogeneous ACDM solution. We consider 
3 models: 

Spherical void, using an LTB model, with the observer 
looking through the centre. The void can be compen- 
sated by a spherical shell of matter, or uncompensated. 
In the compensated case we choose 
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6(r) = 1 -lr- 2 R{r-R)e 3 / 2 - 6 U r - R V R r \R <c r < fi?, 
5 min ~ I -r- 2 (r - 2R) 2 |i? < r ^ 2#, 





and in the uncompensated case 
1 



r > 2R, 



(8) 



r ^ \R, 

\r- 2 R(r - fi) e 3/2-6[(r-fl)/fl] 2 1 R ^ r ^ R 
{ r^R. 

(9) 

(We only consider growing modes, i.e. a uniform big bang 
time.) 



3 



< 



< 



< 




[3], with the same mass distribution M(r) as the com- 
pensated LTB model. In addition there is a dipole- 
like contribution of the form —(S'/S)cos9 where S = 



FIG. 2: Difference in magnitude (top panels) and density con- 
trast (bottom panels) between the background ACDM model 
and the exact void models. We show results for an observer 
looking through spherical compensated (top), quasi-spherical 
(middle, shown for both horizontal - blue, dashed - and ver- 
tical - red, dotted - lines of sight) and cylindrical (bottom) 
voids. The insets show the equatorial density contrast, with 
each square 160 Mpc across. In the spherical case we also show 
the uncompensated case (red dotted line) discussed in Fig. [I] 
The standard weak lensing prediction is shown in green in 
each case, which predicts the wrong sign and amplitude of 
the effect. 



luasi- spherical void, with a mass concentration off to 



for r <2R and S = (2R)-°"e-°" 



one side, with the observer looking either through the 
concentration, or along a line of sight not containing 
the concentration. We use the Szekeres type I model of 



r -0.99 e -0.99T-/(2i?) 

const, for r > 2R. This generates a compensated inho- 
mogeneity extending to r = 2R. 

Cylindrical void with the observer looking along the sym- 
metry axis, using a Szekeres type II model [5]. The den- 
sity profile along the symmetry axis is compensated and 
extends to about r = 2R. 

We standardize the voids in each model radius 50 Mpc, 
with centre at z = 0.1. This corresponds to a slightly 
different comoving distance in each model, and the red- 
shift of the near and far sides are slightly different. The 
depths we choose are d m j n = —0.95 for the spherical and 
quasi-spherical cases, and <5 m i n = —0.8 for the cylindrical 
case. Our model void size is typical in the SDSS survey, 
where void radii are in the range 5 to 135/j^ 1 Mpc [3]. 
Figure [2] shows the change in distance modulus for each 
type of void. In all cases we see the same qualitative be- 
haviour as in the perturbative case: a relative dimming 
of objects on and near the closer side of the void, through 
to a relative brightening of objects located on and near 
the far side of the void. Details of the signal depend on 
the void shape and the nature of compensating regions, 
but the usual weak lensing prediction ^ is always com- 
pletely wrong, unless the source is far from the void (cf 
the insets in Fig. [T]) . 

Thus the predictions of the relativistic perturbative 
magnification persist in the nonlinear regime and are 
generic for different void configurations. But nonlinear 
effects can be large, and the linear relativistic analysis 
can be wrong by more than ~ 20% compared to an exact 
spherical void model. This is shown in Fig. [3j which also 
includes the error for unvirialized over-densities. 

Discussion Our results illustrate a general princi- 
ple: the measured magnitude of astronomical objects de- 
pends not only on internal properties of the source and 
statistics of large-scale structure, but also on environ- 
ment and the nature of inhomogeneity along the line of 
sight, which causes (de-)magnification. (See also [5] and 
[7|.) We have used the full relativistic perturbative anal- 
ysis to show that magnification of objects located in and 
near cosmic voids is dominated by the Doppler term ([7| , 
which overwhelms the lensing effect ^ - and is of op- 
posite sign near the far edge of the void (see Fig. [lj . In 
other words, we have uncovered an anti-lensing effect for 
sources near the far side of voids. The usual weak lens- 
ing analysis fails completely to predict the foreground 
dimming and the background brightening from a void. 
Using exact solutions to model different voids, we have 
shown that this qualitative behaviour persists for non- 
linear voids and for different shapes and ray directions 
(Fig. [2]) . Nonlinear corrections to the relativistic linear 
predictions can be large (Fig. [3]). 

The failure of standard weak lensing for objects in and 
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maximum magnitude difference between fractional difference between the full 

lensing approximation and exact result linear approximation and exact result 
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FIG. 3: (Left) Accuracy of the lensing approximation compared to the full relativistic approximation for the magnification 
near the far edge of a spherical compensated void (S < 0) and a unvirialized over-dense lump (8 > 0). (Right) Comparison of 
the full linear approximation to the exact result: maximum differences of > 20% are seen for deep - but nevertheless realistic - 
voids. Voids/lumps of radius fOO, 50, 25, and 10 Mpc are used. 



near under-dense void regions extends also to over-dense 
lump regions, provided that they are unvirialized. Com- 
mon to both cases is the coherent flow into or out of the 
region, which sources a large velocity contribution. The 
key difference is that the lump occupies a much smaller 
(and shrinking) volume. As the lump continues to con- 
dense it tends to virialize and the correction then dies 
away, whereas the correction for the void increases with 
time. For a source near a cluster of galaxies, there is 
no large net inflow or outflow since the structure is close 
to virialized, and so the usual weak lensing analysis is 
accurate. 

The relativistic linear analysis is accurate for comput- 
ing the velocity contribution due to large-scale structure, 
provided the voids have S m [ n > —0.2 (see Fig. [3]). The 
velocity contribution from large-scale structure was esti- 
mated for the angular power spectrum at fixed redshift 
in [I] - but without taking account of the nonlinear ef- 
fects from voids that we have identified, which we have 
shown amplify the effect. The velocity contribution CJ 
was predicted to exceed the usual Cf for z < 0.2, to be 
about 50% for z sa 0.5 and to be negligible for z > 1. The 
nonlinear corrections that we have identified due to voids 
with S m - ln < —0.2 will introduce a systematic error into 
the perturbative calculation. A key question is: how to 
estimate the nonlinear void correction and thus correct 
for this systematic in galaxy surveys? A similar question 
applies also to supernovae magnitudes [5J. 

While the void effect on convergence is large and of 
the opposite sign to what is expected from lensing, the 
same is not true for lensing shear. The Doppler effect 



will move the redshift of a sheared source, which is a 
correction that should be taken into account in making 
density maps from shear data. However, the sign of the 
shear is unchanged, as the Doppler term does not distort 
an object's shape, but only changes its inferred size. As 
pointed out in [T|, measurement of the convergence and 
lensing shear will provide a powerful probe of peculiar 
velocities. 

Could this effect have already been detected? It is 
reported in [8] that the scatter in the Hubble diagram 
of SNIa magnitudes depends on the environment of the 
host galaxy. They find that galaxies which have a lower 
star formation rate are associated with SNIa which have 
a smaller scatter in the Hubble diagram (at 2 — 3<j). As 
lower star formation rates are associated with galaxies 
in clusters compared to those in the field - the latter 
being more likely to be on the edge of a void - it may 
be that this extra scatter is due to the large relativistic 
anti-lensing effect we have described here. This deserves 
further investigation. 
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